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THE  PRESENT  CALCULATION  METHOD  OP  ADHERING  LAYERS  ON  AIRCRAFT  WINGS 
Bao  Hanling 

With  smaller  and  smaller  tolerances  In  aircraft  design,  there  are  higher 
and  higher  requirements  projected  for  aircraft  performance.  Like  the  wind 
tunnel,  ccnputers  are  also  expected  to  be  an  experimental  tool  In  quanti¬ 
tatively  forecasting  aircraft  performance.  Thus,  it  bee  ones  more  and  more 
urgent  to  calculate  analytically  the  viscous  adhering  layer  and  its  Inter¬ 
ference  effect  on  non-  viscous  flow.  Since  the  1970s,  considerable  efforts 
have  been  given  In  solving  three-dimensional  adhering  layers.  It  is 
considered  that  methods  of  adhering  layers  and  successive  substitution  of 
potential  flow  can  be  used  to  obtain  the  runerical  solution  of  the  adherance 
problem  to  satisfy  engineering  accuracy.  In  the  aeronautical  field,  first 
of  all  the  calculation  of  three-dimensional  aircraft  wings  must  be  solved. 
Since  1975*  almost  all  aircraft  corporations  In  advanced  countries  have 
developed  their  own  computer  programs  for  calculation  of  the  three-dimen¬ 
sional  adhering  layer  of  aircraft  wings.  However,  since  the  engineering 
simulation  of  turbulent  flow  Is  still  not  accurate  enough,  and  there  are 
Insufficient  data  from  wind  tunnel  and  flight  certification,  considerable 
gaps  exist  for  the  application  requirements.  Up  to  now,  development  toward 
perfection  is  still  underway. 

At  present,  considerable  eqphasls  Is  still  given  to  the  development 

of  adhering  layer  calculation  Internationally.  Since  the  first  Turbulent 

Flow  Adhering  Layer  International  Conference  was  convened  at  Stanford 
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University  In  1968,  similar  activities  have  taken  place  several  times.  Every 
few  years,  the  United  States  Air  Force  and  West  Germany  convene  an  adhering 
layer  technical  exchange  conference.  In  Europe,  a  permanent  organization, 
the  Viscosity  Committee,  has  been  founded  to  promote  and  conduct  the  evalua¬ 
tion  of  the  calculation  method  of  shearing  flow. 

As  is  well  known,  the  primary  problem  of  calculating  the  turbulent  flew 
adhering  layer  is  to  make  models  for  turbulent  flow  phenomena.  At  present, 
there  are  two  major  categories  of  methods  used:  field  simulation  and  total 
volume  simulation.  Field  simulation  is  called  the  differential  method,  and 
total  volume  simulation  is  called  the  integration  method.  In  either  method, 
considerable  experience  should  be  related  on  to  completely  solve  the  equations 

Differential  Method 

Generally,  the  differential  method  solves  the  partial  differential 
equations  based  on  the  adhering  layer  equations  including  the  appropriate 
Reynolds  stress  hypothesis.  According  to  different  hypotheses,  the  equations 
can  be  divided  into  two  types:  (1)  one  or  more  differential  equations  to 
directly  represent  the  hyperbola  type  of  stresses,  and  (2)  the  parabola  type 
with  the  use  of  vortex  viscosity  to  represent  stress.  The  degree  of  complex¬ 
ity  of  the  second  type  of  equations  is  determined  by  the  description  modes 
of  characteristic  length  and  characteristic  velocity  of  vortex  motion.  If 
algebraic  means  are  used  to  describe  the  characteristic  length  and  velocity, 
this  is  called  the  model  of  zero  equations.  If  only  the  characteristic 
velocity  is  described  with  differential  equations,  this  is  called  the  single 
equation  model.  If  both  the  characteristic  length  and  the  characteristic 
velocity  are  described  by  differential  equations,  this  is  called  the  double 
equation  model.  Of  course,  it  is  feasible  to  develop  a  multi-equation 
model  but  the  greater  calculation  time  required  cannot  be  compensated  by 
accuracy;  therefore,  it  is  seldom  used.  In  the  method  of  vortex  viscosity, 
the  zero  equation  model  is  extensively  used  at  present,  because  from  the 
study  of  the  two-dimensional  adhering  layer,  accuracy  using  the  single 
equation  model  has  no  significant  improvement,  while  the  double  equation 


model  Is  too  complicated  (in  the  region  nearing  the  separating  flow,  the 
double  equation  model  can  often  give  mere  rational  results).  In  the 
aerodynamic  field,  the  zero  equation  model  with  wide  applications  Includes 
the  Cebec i-Smith  model  and  Mitchell  model.  The  former  model  is  the  dual-layer 
model:  the  danping  factors  of  longitudinal  direction  pressure  gradient  and 
wall  surface  mass  jetting  effect  are  included  in  calculation  of  the  internal 
adhering  layer.  In  the  external  layer,  consideration  is  given  that  the 
vortex  viscosity  is  proportional  to  the  displacement  thickness  of  the  two- 
dimensional  form,  and  also  the  Influence  of  alternating  turbulent  flows. 

The  Mitchell  model  is  a  single  layer  model  composed  of  a  single  hyperbola 
tangent  function  with  multiplication  of  a  simple  Von-Trest  damping  factor. 

Ihe  Cebeci-Staith  model  is  more  rational  because  mare  factors  are  considered; 
however,  this  model  is  more  ccnplicated  than  the  Mitchell  model.  Whatever 
the  type  of  algebraic  model,  it  is  only  determined  by  the  local  gas  stream 
parameter;  i.e.,  the  "upper  stream  historical  effect"  of  the  turbulent  flow 
is  not  considered.  The  upper-stream  effect  can  only  be  revealed  by  vising 
differential  equations  to  describe  the  turbulent  flow  regime. 

Quite  a  few  researchers  often  introduce  the  isotropic  hypothesis  of 
the  vortex  when  using  the  vortex  viscosity  model  in  the  three-dimensional 
state.  In  other  words,  values  of  vortex  viscosities  in  various  directions 
are  the  same.  Actually,  as  early  as  197*1  East’s  experiments  indicated  that 
this  hypothesis  is  actually  without  rational  basis.  Not  long  ago,  papers 
written  by  Tassa  et  al.  revealed  that  the  non-isotropism  of  vortex  viscosity 
has  an  Influence  that  cannot  be  neglected  in  calculating  the  adhering  layers. 

Integration  Method 

Generally,  the  integration  method  Integrates  the  main  equations  in 
advance  along  the  thickness  direction  of  the  adhering  layers.  Thus,  the 
original  parabola  type  equation  with  three  variables  becomes  a  hyperbola 
equation  with  only  two  variables.  After  integration,  the  Reynolds  stress 
term  becomes  the  coefficient  of  friction.  Ihis  equation  is  called  the 
moment  integration  equation  of  the  adhering  layer.  As  proved  by  considerable 
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practice,  a  good  integration  method  should  be  coordinated  with  appropriate 
differential  type  auxiliary  equations.  The  blending  equation  is  a  relatively 
good  auxiliary  equation  currently  used  for  three-dimensional  calculation. 

The  method  of  using  this  equation  is  called  the  blending  method. 

The  rational  behind  the  blending  method  can  be  explained  in  the  follow¬ 
ing:  the  reason  the  energy  of  the  adhering  layer  of  the  turbulent  flow  is 
not  dissipated  is  that  the  energy  formation  process  of  the  turbulent  flow 
(Just  the  opposite  of  the  dissipation  process)  occurs  in  the  adhering  layer; 
this  is  the  so-called  abrupt  phenomenon  of  turbulent  flow.  The  energy 
required  in  the  formation  process  is  provided  by  the  interaction  of  free 
flow  and  turbulent  flow;  i.e. ,  the  fluid  of  the  non-turbulent  flow  (of  the 
free  flow)  at  the  boundary  of  the  adhering  layers  is  continuously  rolled 
into  the  adhering  layers;  this  process  is  called  blending.  The  momentum 
is  transmitted  into  the  fluid  of  the  turbulent  flow  of  the  adhering  layer 
by  collision  between  fluid  masses.  Therefore  the  fluid  mass  entering  the 
adhering  layer  during  the  blending  process  can  more  or  less  reveal  the  total 
volume  characteristic  of  the  turbulent  flow  phenomenon  in  the  adhering  layer. 

The  blending  coefficient  Cg  is  the  most  inport  ant  parameter  in  the 

blending  method;  the  coefficient  represents  the  dimensionless  flow  rate  of 

the  net  input  into  a  controlled  cross-sectional  area  of  the  adhering  layer. 

For  different  derivation  methods  of  Cg,  there  are  three  branches  of  the 

blending  method:  1.  the  Smith  method:  C  is  the  algebraic  form  obtained 

E 

empirically;  2.  The  Cousteix  method:  CL  is  the  algebraic  fora  obtained 

through  similar  solutions;  and  3*  The  Stock  method:  C  is  derived  by  an 

E 

enpirical  differential  equation  (thus,  consideration  is  given  to  the 
’’upper-stream  historical  effect”  of  the  blending  process). 

In  the  integration  method,  only  the  total  volume  parameters  of  the 
adhering  layer  are  given,  such  as  the  coefficient  of  friction  at  the  wall 
surface.  The  flow  details  of  the  flow  field  are  not  provided.  Since  the 
fundamental  basis  of  the  integration  method  is  not  as  sound  as  it  should  be, 
the  integration  method  is  inferior  to  the  differential  method  in  adapting 
to  different  types  of  flew.  At  the  present  stage,  however,  the  integration 


method  is  as  good  as  the  differential  method  with  fast  calculation  speed, 
quite  high  cost  effectiveness,  and  is  more  adaptable  for  application  in 
calculations  requiring  repeated  and  successive  substitution  for  the  adhering 
layer  and  the  potential  flow.  Therefore,  in  many  countries  research  and 
development  are  still  stressed  on  this  type  of  method. 

Table  1  lists  methods  and  types  for  calculating  the  three-dimensional 
adhering  layers  of  an  aircraft  wing. 

Table  1. 


in 

*  (b)  « 

IQCTW 

iaanEgi 

n 

A.Tatu9  (e) 

«*«*«»«* (f) 

AIAA-U-0224  | 

H 

J.F.N**k*(e) 

«»£**«;**  (K*)(g 

■ZSE2EH3SHBI 

ft 

J.F.N«k*(e) 

wwetjtmim 

Lockheed  Rep.  LGjfERoiM 

4 

T.Ccbtcitp  (g) 

Mtmttnmma  (f ) 

NASA  CR-2nr « (e) 

§ 

I.D.McUu 

- (e~) 

**«««»««  (f) 

AIAA77-S 

4 

I.F.F.LUdfcontft 

••e  Msmis* 

NLR  TRramU  V  (e) 

t 

J.ComMx* 

ONERA  TR  N0it75-4S 9(e) 

$ 

W.KorSalla 

**«*«»«*  (f)- 

AIAA -77-30*9  (e) 

D 

A*K*  RmI4|iF 

S*tt*K#««a(l) 

A1AA  I. Volts  Not 

D 

»aw*(k) 

ARC  R&.M  arts 

Eszzam 

»M*(lc) 

Doraitr  GMBHn/llB  9  (e) 

L-J 

D.B.Speldiaf 

Limn  EEEBfl 

Coa peter  llethod  |a  Applied 
Mechanics  end  EaciaeeriatVolS.Nl 

Key:  (a)  Sequential  number;  (b)  Authors;  (c)  *fype  of  methods; 

(d)  Source  of  papers;  (e)  et  al.;  (f)  Algebraic  model  of  vortex 
viscosity;  (g)  Reynolds  stress  equation  method  (exposed  type); 

(h)  Reynolds  stress  equation  method  (concealed  type);  (i)  Double 
equation  model  of  vortex  viscosity  (incompressible);  (k)  Blending 
method;  (1)  Double  equation  model  of  vortex  viscosity  ;  (j)  Reynolds 
stress  equation  method  (incompressible). 

Coordinate  System 


One  of  the  important  techniques  in  calculating  three-dimensional 
adhering  layers  is  the  selection  of  a  coordinate  system.  A  good  coordinate 


system  should  adapt  to  a  random  plane  shape  of  an  aircraft  wing  and  simplify 
calculation  as  much  as  possible.  Generally,  these  two  requirements  are 
contradictory.  In  the  rectangular  coordinate  system,  the  expression  mode 
of  equations  is  sinple  but  Inconvenient  in  usage.  Most  researchers  using 
the  differential  method  are  accustomed  to  the  non-orthogonal  curvilinear 
coordinate  system  adhering  onto  a  surface  of  substance:  the  surface  of  an 
aircraft  wing  is  divided  into  lattices  by  percentage  ratios  of  chord  length 
and  span  length  into  coordinates  of  x  and  z  directions;  the  direction 
perpendicular  to  the  surface  of  the  aircraft  wing  is  the  y  coordinate  with 
standardization  treatment  in  order  to  reduce  variations  on  the  thickness 
of  the  adhering  layer  and  variation  of  flow  variables  between  two  adjacent 
wing  cross-sections,  so  that  the  step  length  of  the  difference  can  be 
relatively  long.  Hie  integration  method  used  abroad  adopts  a  random  curve 
coordinate  system  without  restricting  the  plane  shape  of  the  aircraft  wing, 
but  the  method  has  very  long  calculation  formulas.  In  China,  a  blending 
method  program  adaptable  to  calculation  of  the  adhering  layer  of  a  random 
plane  shape  has  been  used.  In  this  method,  equations  are  written  in  the 
rectangular  external  streamline  coordinate  system.  When  equations  are  being 
solved,  the  equations  are  converted  into  the  equal  percentage  ratio  chord 
coordinates;  thus,  the  method  becomes  relatively  sinple. 


Difference  Format 

In  both  the  differential  and  integration  method,  in  the  last  stage 
difference  method  solving  for  differential  equations  is  required.  In  solving 


parabola  type  equations  in  the  differential  method,  customarily  a  three- 

point  or  two-point  (with  second  stage  accuracy)  box  type  concealed  type 

difference  format  is  adopted  along  the  directions  of  the  chord  and  thickness; 

in  the  span  direction,  the  exposed  type  difference  format  is  adopted.  Then, 

in  selecting  difference  elements,  the  concept  of  the  related  region  should 

be  considered.  In  the  region  where  the  span  direction  flow  has  its  direction 

changed,  in  some  methods  there  is  only  first  stage  accuracy  along  the  span 

direction.  In  some  other  methods,  along  the  thickness  direction  of  the 

adhering  layer  and  the  span  direction,  there  are  alternately  adopted  the 

exposed  and  concealed  types;  this  is  called  the  ADI  method.  In  this  method, 
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a  faster  convergence  is  obtained.  For  example,  on  a  V.AX1 1/780  ccnputer,  the 

_2 

average  calculation  time  of  each  lattice  point  is  only  6.6x10  second  in  CPU. 

Both  the  Reynolds  stress  equation  and  blending  methods  are  of  the 
hyperbola  type.  When  the  equations  are  being  solved,  both  the  exposed  and 
the  concealed  types  can  be  used.  Since  a  long  time  ago,  it  was  considered 
that  the  efficiency  of  the  exposed  type  is  relatively  low;  however,  as 
revealed  in  calculations  recently,  the  overall  cost  effectiveness  of  the 
exposed  type  method  is  better  than  the  concealed  type  method  because  it 
avoids  the  time-consuming  successive  substitution  process  though  with  a 
relatively  short  step  length. 

Ccnparison  of  Methods 

Although  many  calculation  methods  relating  to  the  adhering  layer  of  an 
aircraft  wing  were  developed  over  many  years,  up  to  now  we  still  cannot  say 
which  calculation  method  and  code  have  sufficient  accuracy  through  sufficient 
operation  conditions  of  flew  along  aircraft  wings.  It  ».  13  published  at  AIAA 
(1980),  the  results  of  four  calculation  methods  for  three  different  plane 
shapes  of  aircraft  wings  (2,  3,  4,  and  5  in  Table  1)  used  by  three  major 
aircraft  corporations  in  the  United  States  (Lockheed,  Douglas  and  Boeing). 
Calculations  were  conducted  on  CDC  7600  computers.  Refer  to  Table  2  for 
the  calculation  time  used  on  the  number  of  dividing  lattices  and  the  lattice 
nodal  points  on  the  surface  of  each  aircraft  wing. 


As  revealed  by  calculation  results,  there  are  considerable  differences 
between  various  methods,  especially  between  the  exposed  and  concealed  type 
methods.  This  reveals  that  the  solution  of  adhering  layer  equations  not 
only  is  very  sensitive  to  different  mechanical  models,  but  also  to  different 
mathematical  processing  methods.  Hie  conclusion  of  the  report  is  that  it 
is  urgently  required  to  have  reliable  experimental  results  to  further 
clarify  which  methods  are  the  better  ones.  Besides,  the  report  also 
considers  that  for  aircraft  wings  of  some  plane  shapes,  sometimes  the 
difficulty  of  non-cc  merger  in  successive  substitution  may  occur  in  the 
calculation  process  of  mo  algebraic  model  methods. 


Table  2. 
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Total  number  of  nodal  points  *  span  direction  x  chord 
direction  x  normal  direction. 

Nunber  of  nodal  points  (at  the  surface  lattice  of  the 
aircraft  wing)  N  *  span  direction  x  chord  direction. 

Total  CPU  time  =  tx^N  . 
y  P 

Key:  (a)  Code;  (b)  Aircraft  wing;  (c)  Second;  (d)  Nash 
exposed  type;  (e)  Nash  concealed  type;  (f)  McLean;  (g) 
Cebeci. 


In  1981,  a  shearing  flow  calculation  method  evaluation  team  of  the 
European  Viscous  Flow  Conmittee  used  a  35°  sweptback  aircraft  wing  and  flow 
M  number  0.5  in  a  calculation  example  to  conduct  unified  evaluation  on  eight 
established  calculation  methods  for  the  adhering  layer  of  turbulent  flow 
In  European  countries.  Seven  of  the  eight  calculation  methods  used  methods 
5  through  11  given  in  Table  1;  another  method  is  the  Brest  method  of  two- 
dimensional  stress  equations.  The  evaluation  results  of  the  calculation 
methods  are  similar  to  the  aforementioned  conclusions,  in  that  the  differen¬ 
tial  method  is  no  better  than  the  integration  method  with  higher  accuracy. 
Therefore,  considerable  work  should  be  conducted  to  reach  the  practical 
degree  far  calculating  the  three-dimensional  adhering  layers  of  aircraft 
wings. 


In  1982,  Tassa  et  al.  in  the  United  States  published  a  relatively 
encouraging  calculation  result;  his  method  is  similar  to  the  Cebeci  method 
except  the  calculation  on  the  influence  of  isotropism  in  vortex  viscosity 


and  adoption  of  the  ADI  technique.  In  other  words,  the  Tassa  method  utilizes 
all  the  good  points  of  the  currently  available  differential  method.  As 
revealed  in  results,  the  calculation  moment  thickness  and  chord  direction 
displacement  thickness  match  the  experimental  results  quite  well;  however, 
the  result  of  the  span  direction  displacement  thickness  is  not  as  good  and 
the  calculation  examples  are  only  limited  to  the  flow  field  of  the  central 
region  of  the  half  wing-span.  Further  checking  is  required  on  the  accuracy 
of  the  entire  aircraft  wing. 

Several  Viewpoints 

1.  In  the  calculation  result  affecting  the  three-dimensional  adhering 
layer  (especially  the  result  of  lateral-direction  parameters),  one  of  the 
most  outstanding  problems  in  accuracy  is  the  realization  of  the  direction 
field  of  Reynolds  stress  or  the  non-isotropism  of  turbulent  flow.  At  present, 
in  some  methods  the  non-isotropic  hypothesis  used  is  still  far  from  perfect, 
as  it  is  without  theoretical  basis. 

2.  The  integration  method  has  an  advantageous  aspect  because  it  only 
solves  for  the  total  volume  and  avoids  the  difficult  problem  of  the  stress 
direction  field.  However,  it  is  difficult  to  process  the  lateral-direction 
velocity  cross- sect ional’~st ate ,  and  similarly  there  is  considerable  difficulty 
in  accurate  simulation  of  lateral-direction  parameters.  There  will  be 
greater  difficulties  when  reversal  of  lateral  flow  occurs  in  some  regions 

of  the  aircraft  wing. 

3.  When  the  flow  is  nearing  separation,  the  calculation  becomes  very 
difficult  due  to  the  singularity  of  the  adhering  layer  equation  itself. 
Therefore,  frequently  a  cutoff  calculation  is  required  at  a  certain  position 
before  separation.  This  requires  a  criterion  of  determining  three-dimensional 
separation.  This  is  also  one  of  the  important  problems  worth  stressing  in 
present  and  future  studies. 

Jj.  Since  the  input  data  are  quite  sensitive  in  calculating  the  corre¬ 
sponding  potential  flow  of  the  adhering  layer,  it  is  necessary  to  develop  a 


high  accuracy  calculation  method  in  developing  the  potential  flow  for  mutual 
successive  substitution  of  adhering  layer  flow  and  potential  flow. 

5.  The  greatest  obstacle  hampering  the  theoretical  development  of  a 
three-dimensional  adhering  layer  is  the  lack  of  experimental  data  used  in 
certification.  Since  the  interference  effects  of  the  wind  tunnel  wall, 
fulcrum  and  probes  still  cannot  be  sufficiently  eliminated,  at  present 
considerable  technical  difficulties  still  exist  in  measuring  the  adhering 
layers,  especially  in  a  high  air  speed  wind  tunnel. 

i 

6.  Studies  in  the  1970s  on  the  large  scale  vortex  phenomena  and  the 
abruptly  occurring  phenomenon  of  turbulent  flow  at  the  laminar  bottom  layer 
of  the  adhering  layers  can  be  said  to  be  a  historical  discovery  in  studies 
of  turbulent  flow.  In  the  1970s,  the  so-called  sub-scale  lattice  method 
appeared;  in  the  method,  the  three-dimensional  scale  used  is  considerably 
smaller  than  the  largest  vortex  scale,  but  the  three-dimensional  scale  is 
still  considerably  greater  than  the  smallest  vortex  scale.  Therefore,  the 
influence  of  large  scale  vortexes  can  be  revealed.  Although  the  sub-scale 
lattice  method  is  still  at  the  initial  stage,  this  method  shows  a  new 
direction  in  studies  of  turbulent  flow. 

7.  Before  the  1970s,  the  calculation  of  three-dimensional  adhering 
layers  of  an  aircraft  wing  was  almost  a  blank.  In  this  short  span  of  ten 
years,  more  than  10  methods  have  appeared.  Although  these  methods  are  not 
entirely  satisfactory,  it  can  be  said  that  the  mathematical  difficulties 

of  solving  the  equations  are  nearing  a  complete  solution.  Only  in  mechanics 
are  there  still  many  unknown  factors  waiting  for  researchers’  creative 
explorations  in  continued  studies. 


